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-Tunnel Test: Problem in Superae 


Hsue-shen Tsien 
Guggenheim Aeroanutical Laboratory 


Massachusetts Institute of Technology 


Wind-tunnels are, perhaps, the most useful tool in aerodynamic in- 
vestigations and certainly have contributed much in the modern develop- 
ment of fluid mechanics, It is thus natural, when one turns to a new field 
of aerodynamics, the aerodynamics of rarefied gases or Buperaerodynamics, 
that one should think of using the wind tunnel again. Only here it has to 
be adopted to entirely new circumstances and many new problems, both in its 
design and in its operation, appear. It is the purpose of this paper to dis- 
cuss Some of these problems, so as to gain an orientation in the new field 


of experimentation. 


l. Tunnel Design 
To test models in the wind tunnel at its test section, it is of 


primary importance to obtain a uniform stream at the desired temperature, 
pressure and velocity. For subsonic wind tunnels with ordinary pressures, 
this can be achieved without much difficulty. For supersonic wind tunnels 
&t ordinary pressure, the expansion part of the tunnel before the test sec- 
tion or the nozzle is first designed to obtain a uniform stream at its exit 
without considering the effects of the viscosity of air. Then the boundary 
layer along the wall of the nozsle is calculated with the pressure gradient 
thus determined. Finally the displacement thickness of the boundary layer, 


or the needed space for the boundary layer flow at lower velocity is provid- 


ed by making the (0) i determined by the 


calculated amount. (Ref. 1) This design procedure is found to give satis- 
factory nossles for supersonic wind tunnels. 

However when one tries to use the same design procedure to the super- 
aerodynamics wing tunnel, one is immediately confronted with the difficulty 
of extremely large viscous effect, In other words, the boundary layer will 
be so thick as to occupy the main portion of the noszle passage. To demon- 
Btrate this effect, let the length of the test section be L and the width of 
the square test section be b. Then the Reynolds number based upon the oondi- 
tions in the test section ie Re = ኞ where U is the velocity in the test sec- 
tion. If, as a rough estimate, we take the thickness of the boundary layer 
to be zero at the beginning of the test section and equal to 5 value calcu- 
lated by the well-known Blasius formula for a flat plate at the end of the 
test section, then 

8 = 3.65 L/ e (1) 


Now if this boundary layer actually occupies half the tunnel width b/2, then 


b 
$=2 3.63 LNS (2) 


On the other hand, the ratio of the mean free path E and the boundary layer 


thickness $ is known (Ref. 2) to be equal to 


L a 12557 M (3) 


j 2.65 NRe 


where 1 is the ratio of specific heats and can be taken as 1.4, M| is the Mach 


number in the test section. By combining (2) and (3), we have 


Fal ቴ) = 0.0557 M (ፊ) 


D AAA 


m" 


This relation 46 a OJ = a d #8 equal to 2, 


and L/b = 2, the boundary layer will completely fill up the test section, 

if the mean free path is equal to 5.6% of the boundary layer thickness or 
2.8% of the tunnel width. This means that the extremely strong viscous ef- 
fect at low derivities makes the ordinary concept of designing a wind tunnel 
totally inapplicable. 

The extremely thick boundary layer where the velocity increases from a 
small value near the wall to some supersonic velocity at the center of the 
nozzle, also gives subsonic velocities in a rather large portion of the nos- 
zle. Since pressure disturbances downstream can be transmitted upstream in 
subsonic flows, the flow in the test section of a low pressure tunnel will 
be sensitive to changes in the diffuse even if the main stream velocity at 
the center of the nozsle is supersonic. This is of course a new phenomenon 
is superacrodynamics tunnel not found in conventional supersonic wind tunnels. 

The large viscous effects can be also demonstrated by calculating the 
ratio of frictional loss on the walls of the test section and the shock loss 


in the diffuser after the test section. Consider the diffuser to be a straight 


tube of approximately the same cross-sectional area as the test section, then 


the pressure loss due to friction ak, , 28 


AP, = Frictional Force 


y? 
EK ie d 
DATE 
1,328 


Taking (, to be Blasius value or hr 


we have 
1328 (5) 
N Re 


i è Now if is the QJ E Bec then the pressure by 


ideal isentropic compression in the diffuser is vds 
1) the actual pressure rise in the díffuser 1s estimated as that due to a nor- 
mel shock without further recovery, then the actuel pressure rise is 


Therefore the pressure los8 due to shock is 


(6) 


By oombining (5) and (6), the ratio of these two pressure losses is 
| (7) 
Introducing the mean free path ratio given by (3), we rave 
(8) 


| This relation is plotted in Fig. 2. Therefore if the Mach number 1s 2, 
and L/b = 2 as before, then when the ratio ( ) ís 0.056, the ratio of 
frictional loss to shook loss is 0.628. Hence the frictional loss and the 
shock loss is of the same order of magnitude. 
These large viscous effects are fully confirmed by the recent tests on 
the 1" x 1" low pressure wind tunnel at the University of California," The 
test nossle (Fig. 3) was designed for Mach number without considering 


E the viscous effect of the medium. During test, the static pressure on the 


* Experimental work done under contact with the Office of Naval Research. 
The author ia deeply indebted to Professors R. G. Folsom and E. D. Kane 


for permission to use their unpublished results. 
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wall at the exit (COLE D 48 equal to 175 


microns" for the two tests preserited in Figs. 4 and 5, The apparant Mach 
number "M" is the Mach number calculated from the dynamic pressure measured 
by a Pitot tube by using the Rayleigh formula. Since there is the compli- 
cation of large viscous effect in the Pitot tube reading as shown in the fol- 
lewing section, this apparant Mach number is only quantitative and cannot be 
taken as the exact value. However it is apparant from Figs. 4 and 5 that the 
boundary lsyer in the test section is indeed very thick, and fille up the 
whole spice, This large boundary layer thickness mikes this spacco available 
for th: expansion of the central potential flow if it exists, very small. 
Therefore tha maximum Mach number reached at the center of the nozzle is very 
much smaller than the design Mach number of 4. At the lower pressure, the in- 
fluence of slip at the wall is also evident. This has the tendency to make 
the flow more uniform. However the very low Mach number at the test section 
indicates again the strong viscous effeot in converting much of the pressure 
energy into hent energy. 

These elementary calculations and preliminary test results makes it clear 
that for the design of the nozsle end test section for « superaerodynamics wind 
tunnel, it is no longer possible to separate the compressibility effects and 
the viscous effect. In fact, the concept of boundary layer is also of doubtful 
value due to the extremely small Reynolds number encountered, Therefore to aot- 
ually design such a nozzle to obtain the nearest approximation to the ideal uni- 
form flow it will be necessary to use the exact Navier-Stokes equations instead 
of the approximate boundary layer equations. Of course, it may be argued that 
for superaerodynamics, the Navier-Stokes equations for no more exact and addi- 


tional corrections must be added (Ref. 2). However, recent investigations by 


6 
* 1000 microns = 1 mm Hg, one atmosphere = 0,760 x 10 n. 


dde 


R. Schamberg E ©) IE" corrections are small 


in case of slip-flows concerned here and will not essentially alter the flow 
pattern. Hence for 8 first approximation just like the non-viscous isentro- 


pic flow as a first approximation for ordinary supersonic noszles, we can use 


the Navier-Stokes equations. The simplest case to be considered is certainly 
the axially symmetric nozsles. If x is the coordinate in the axial direction, 

the coordinate in the radial direction and u and v are the corresponding 
velocity components, the equations are: 


= dessipation function 


= stresses tensor 


(9), (10), (11) and (12) together with the equation of states 


“RT (13) 
then determines the five unknown u, v, , and T. Of course, the actual process 
of making this calculation will be extremely tedious and some approximation 
method of solution may have to be developed. One possibility would be to adopt 


D) N 
the NP case: We integrate i 


the differential equations once with respect to r and thus only try to sat- i 


È 
isfy the equations "on the average" over the cross-section of the nozzle, ፍፄ 
The "distribution" of u, v over the cross-section will then be set in the dy 
form of a polynonic in r. Initial study in this procedure is already made =» 

-T 


by 8. A. Schaaf (Ref. 4) at the suggestion of the author. 

For ordinary supersonic diffuser, high efficiency of pressure recovery | i 
can be generally achieved by using a long diffuser. However, for superaero- | 
dynamios wind tunnel, due to the extremely large loss through friction, long 
diffusers are undesirable. In fact, the pressure loss can be reduced by us- 
ing & shortest possible diffuser. 


2. Flow Measurement 
The quantities which determine the fjow field are two out of 
the three veríables  , „ T and the velocity components. The quantities 
„ 7, T are related by the equation of states and therefore only two is nec- 
essary for the determination of all three, Generally for wind tunnel work, 
the quantities actually measured are „ „ and  , the magnitude of the velo- 
city. 

For the measurement of pressure, a manometer is used, For ordinary pres- 
Bure, one uses a fluid manometer filled with water, alcohol or mercury. How- 
ever for the extremely low pressure encountered in the superaerodynamic flow, 
some other form of manometer is necessary. One of the most successful type is 
the Pirani gage. The conventional form of Pirani gages has a pressure sensiti- 
vity of about 10? micron. It utilises the change of temperature of a wire 
heated with constant energy caused by a change in the pressure of the gas sur- 
rounding it. The temperature change is measured by the change in the resistance 


of the wire, The wire is looated in & small chamber which is connected to the 
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point of Eg ፈሕ Ed if static pressure 


is to be measured, The question of best design of the connecting tube for 
quick response is studies by 8. A. Schaaf. (Ref. 5). 

To measure the density , the conventional method utilises the differ- 
ence in the velocity of the light rays in mediums of different density. With 
different optical arrangements, we have the shadowgraph method, the schlieren 
method and the interferometer method. However, if the density of the medium 
is very low as the case of superaerodynamic flow, the sensitivity of these 
methods become extremely poor. For instance in case of schlieren method, the 
percentage change in illumination I by passing through a region of thickness 
b is given by 


(14) 


where is the air density at 32°F and 1 atmosphere pressure, and 

is the density gradient normal to the light ray. and are the local length 
and the normal unobsoured width of the light source image perpendicular to the 
knife edge. k is factor of order 1, determined by the particular optical path 
used, Therefore the sensitivity of the schlieren method decreases with the fac- 
tor ( ). Some improvement can be made by altering the quantities f and , 
but practical limitations and diffraction difficulties do not allow the increase 
of sensitivity to satisfactory values. 

A new approach to this problem density measurement is the method of absorp- 
tion. It is found for instance, that oxygen at low pressures shows a strong ab- 
Borption band at wave lengths around 1470? A or ultra-violet light. The per- 
centage absorption is proportional to the number of molecules that meet the light 
ray and is, therefore, proportional to the density of the gas. The measurement 
48 then similar to that of the interferometer method where the density is deter- 


mined. A similar method is the utilisation of the after-flow of nitrogen. These 


The conventional method for the measurement of velocity is through the 


use of dynamic pressure rise in a Pitot-tube. ል straight forward application 


of this method is, however, difficult for rarefied gases. The formula used 
is, however, based upon the neglection of viscosity effects. But for rarefied 
gases the viscosity effect is af great importance as pointed out in the previ- 
ous section, Then the dynamic pressure would be quite different than that 
given by the usual formula. To estimate this effect, lot us consider the case 
of low Mach number so that compressibility effects can be neglected. Then as 
a first approximation, take the flow field around the Pitot-tube as that of a 
source of strength in non-viscous flow of uniform velocity U. The "radius" 


of the tube ል is 
as 


and the stagnation point is located at 


(35) 


The velocity introduced by the source 48 then 


By calculating the viscous stress from this approximate disturbance velocity, 
we have for flow along the axis 


(16) 


Hence ፈያ is the stagnation pressure and the Static pressure, 


(37) 


For rarefied gase 6 ©) Kë of the Reynold's 


number of the Pitot-tube could be of the order of unity. Then the dynamic 


pressure rise -~ is not the usual values $ but a value much less than 
that. In fact previous investigation by Baker (Ref, 7) and F. Homann (Ref. 8) 
indicate that the Reynolds number where a! is the radius of the mouth of 
tube, must exceed 30 in order to reaoh the usual dynamic pressure rise è 

When the velocity of flow is high, we have the added complication due to 
the shock. The conventional Rayleigh formula for Pitot tubes is supersonic 
flow is based upon the assumption of very thin shock wave ahead of the Pitot 
tube, Now the thickness of the shock is proportionnl to the mean free path. 
Hence in rarefied flows, the thickness of the shock will be so. increased as to 
cause interference with flow in the neighborhood of the Pitot-tube. This to- 
gether with the viscous effect mentioned in previous paragraph definitely show 
the inapplicability of the Rayleigh formula for supersonic velocity of rarefied 


gases. 


3. Hot-Nire Anemometer 

With the great complications in applying the conventional velo- 
city measuring device to superaerodynamic flows, one is naturally led to the 
thought of other avenues of approach. One possibility is the use of hot-wire, 
If the wire diameter is of the order of 0.0001 inches, and if the pressure of 
the gas stream is approximately 100 microns, the ratio of the mean free path 
to the wire diameter will ^e approximately 180. Therefore the flow around the 
wire is definitely the free molecule flow (Ref. 2). ሸ6 have thus a simple phy- 
sical situation, which is an improvement over the rather uncertain circumstances 
of mixed dynamic and viscous effects for the measurement of velocity by Pitot 
tube, It thus seems worthwhile to explore this possibility by a trial calcula- 


tion of the performance of such a hot-wire anemometer. 


HT AE MR 


ይ macroscopic velocity ሀ and a Maxwellian molecular velocity distribution, 


the translational energy of molecules incident upon the unit area is, 


(18) 
where c? = 2 RT, T temperature of the gaa stream and erf is the error function, 
Now let r be the radius of the hot-wire. Then the total energy E, incident up- 
on a unit length of the wire is the sum of translational energy «nd internal 
energy. If 0, is the specific heat at constant volume, this total energy per 
unit length of wire is 


(19) 
The integrals in equation (19) can be expressed in terms of tabulated functions, 
(see Appendix) Thus 


(20) 


where (21) 


(22) 


The Io and I, are the modified Bessel functions of the first kind of orders sero 
and one respectively. The functions F} and Fg are tabulated in Table 1. 
If T, 48 the wall temperature, and x the accomodation coefficient, the 


naf A ` 


difference NON t upon AR &nd the energy 


E, carried by the molecules re-emitted from the surface is givon by 


ጻ-ጁ" (- * 


whore E, is the energy that would be carried away by the molecules if the 
re-emission were at ths temperature T, of the wire. Therefore 


ጂ-ጁ። cr ኔሮ + (eR + 0 ( 20 U (7) 


+ 1 . 0 0-1) (23) 


This difference of energy is then the net energy input to the wire per unit 
E. length of the wire by the air stream. ^u 


If 1 18 the electric current heating the wire and is the resistance 


of the wire per unit length at the wire temperature, the heat input per unit (o 
length of wire by the heating current is y » Heat is lost from the wire L 7 
by radiation. If 13 Stefan-Boltamann constant and is the emissivity of d | 
the wire surface, the radiation hent loss per unit length is e us int 
Therefore if the wire has reached a steady condition, the heat balance re- 9 
quires d 
4 [ 
T # 
VI déi 
4 (24) We 
| AR 
if This equation can be put into somewhat simpler form by using the relation that DL 


R= Cp = Gel -1) (25) 


45 „ and the correspondin, temperature coefficient of the re- 
sistance is men the resistance can be expressed as 


6 1+ (ET) 


Then trom equation (26) 


L TP iii (28) 
T 1 
a o 
Now introduce ነ, as the reference density and i, as the reference heating 


current, then equation (24) can be written as 


(29) 


The partiouler values of the reference temperature To, the reference den- 
sity o and the reference current i, are not yet fixed. We fix these quanti- 
ties now by requiring that 


1,51 (30) 


OPT 


(23) 


This is then the performance equation of the hot-wire in free molecular flow. 
Now let us investigate in greater detail the case of a bright platinum 
wire. To satisfy equation (30), 


To = 492° R. 
The value for and can be taken to be 0.08 and 0.90 respectively. Then 
equation (31) gives the corresponding pressure P, for , and T, as 


Let the radius r of the wire be 0.0001 inch. Then equation (32) gives the 


reference heating current i, as 


io = = 0.274 milliampere. 
where the resistivity of the platinum is taken as 10.96 x 10$ ohms-cns, There- 


fore the order of magnitude of the different quantities is entirely satisfactory. 
If the wire is used with a constant heating current, then equation (33) can 
be used to calculate the relation between the resistance ratio and the velocity 
ratio ( ፪ ) at constant air stream density and temperature. This is done for 
La =2 mad Va si and the result is given in Fig. 6. It is seen 
that the sensitivity of the instrument is good. Of course, the behavior of the 
hot-wire anemometer will be actually determined by calibration for any experiment. 


„ The author is indebted to Mr. L. Mack for the numerical computations. 
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for the model test to be similar to the prototype, the model must be made of 


same surface material as the prototype, and the fluid must be the same, and 
furthermore the following parameters must be the same: 

(1) Reynolds number Re 

(2) Mach number M? 

(3) free stream temperatune T° 

The radiation heat loss from the surface 38 equal to per unit ares. 

However, if the model is surrounded by the walls of the test chamber, there is 
also an heat input due to radiation from walls of the test section to the model, 
Let us call thís quantity + Then the net heat loss per unit area of the sur- 
face of the model is . This quantity can be rendered non-dimensional 
by dividing it by » Call this new parameter , then 


(34) 


For the prototype, the heat from the walls of the test chamber is absent but K 
there may be solar radiation and the radiation from the earth and surrounding 
atmosphere, (Ref. 0). Denoting this amount by , then the parameter for 


the prototype is 
(35) 
In order for the flow to be same also with respect to the radiation heat transfer, 
(36) 4 


Because of the previous conditions on the Reynolds number and free stream tem- 


perature, (36) is the same as 


(37) 


Since the Nd i EE. K by the accommodation 


coefficient as shown by equation (29), it will be necessary to find mater- 


4818 which cen hold this coefficient constant for a considerable period of time 
So that no frequent calibration is required. However the present analysis seems 
to indicate the feasibility of such an instrument for measurements in rarefied 
gases and further research is definitely desirable. 


4. Parameters ot Flow 
The two parameters that are directly connected with the flow field 
are the Reynold's number Re, defined as 


ጩ 


Res”, 


where is the kinematio viscosity, L the typical linear dimension of the body; 
and the Mach number MO of the free stream. This is true even for slip flows and 
free molecule flows due to the fact that the ratio of mean free path to the typi- 
cal dimension can be also expressed in terms of the Reynold's number and the Mach 
number. 

However, as the pressure or density is reduced, the solid boundary of the 
flow enters actively into the flow conditions by requiring not only that the 
microscopic stream velocity be tangential to the surface but that the interaction 
of the molecules and the wall be considered and that the radiation of energy to 
and from the wall be taken into account. The interaction of the molecules with 
the wall is so far expressed through the fraotion of molecules that are dif- 
fusely re-emitted from the wall, and the accommodation coefficient . Itis 
known that both and are funotions of the temperature of the wall and we have 
reasons to believe that they are also functions of the molecular velocity distri- 
bution. Therefore the interaction of the molecules with the wall ia the same 


only if the wall temperature, the gas temperature and the Mach number of the gas 


0/0 ፻; 72 ያ È 


i 0 1.77245 9 di 
i 0.2 1.73751 0.07020 ' 

0.4 1.63880 0.27269 

0.6 1.49248 0.58560 

0.8 1.32021 0.97843 

1.0 1.14328 1.42053 

1.2 0.97825 1.88555 

dek 0.83480 2.35492 

1.6 0.71628 2.81812 

1.8 9.62252 3.2717 

2.0 0.54683 3.71356 

dua 0.48790 4.14672 

2.4 0.44090 4.57300 

2.6 0.40268 4.99395 

2.8 0.37100 5.41117 


3.0 0.34420 5.87498 


Evaluation of the Functions 71 and 72 


For the function Fi, 


whore Io is the modified Bessel function of first kind and under sero. The 


last step is made possible by the substitution = 
For the function Fa, 


By the definition of the error function, we have 


This form can be simplified by partial integration. Thus 


linear dimension of the prototype. This means that tie wall temperature of 
tha test chamber must be 86 controlled that satisfies (37). 

This set of rather strict similarly rules for model testing in super- 
aerodynamic flows is certainiy difficult to satisfy. In what way the rules 


can be relaxed is the problem oi future research. 
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where I) is the modified Bessel function of firat kind and order one. 


June 12, 1948 


Mr. Robert R. Dexter, Secretary 
Institute of the Aeronautical Sciences 
2 East G4th Street 

New York 21, N.Y. 


Dear Mr. Dexter: 


Enclosed are two manuscripts which I am submitting 
for consideration of publication in the Journal of the 
Aeronautical Sciences. 


The first paper, "Wind-Punnel Testing Problem in 
Supsraerodynamics", is the revised form of my paper presented 
in the Annual Meeting. I apologize for the long delay in getting 
the manuscript into this form. 


The second paper, "Two Dimensional Airfoile in 
Hypersonic Flows", is written by Mr. Richerd D. Linnell and 
1s teken from his Thesis at the Massachusetts Institute of 
Teohnology. 


Very sincerely yours 


H. S. Tsien 
Professor of Aerodynamics 
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2 EAST 64TH STREET 
NEW YORK 21. N. Y. 


July 20, 1948 


Dr. H. $. Tsien 

Guggenheim Aeronautical Leboratory- 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


Dear Dr. Tsien: 


We are now preparing your paper, "“ind-Tunnel 
Testing Problems in Superaerodynamics," for publica- 
tion, and there does not seem to be an abstract or 
summary included. Since we feel that some sort of 
an abstract should be published with the paper, I 
should appreciate your sending one as soon as possible. 


Very truly yours, 


Tomis N re 


Berneice H. Jarck 
Associate Editor 


OFFICIAL PUBLICATIONS 
AERONAUTICAL ENGINEERING REVIEW JOURNAL OF THE AERONAUTICAL SCIENCES AERONAUTICAL ENGINEERING CATALOG 


| July 25, 1958 ብ 


Miss Berneice H. Jarok 

Associate Editor 

Journal of the Aeronautical Sciences 
2 East 64th Street 

New York 21, N.Y. 


Dear Mise Jarck: 
Enolosed ís the Summary for my paper, "Wind-Tunnel 


Testing Problems in Superaerodynamios". This is to answer 
your request in your letter of July 20, 1948. 


Very sincerely yours 


He 8. Tsien 
Professor of Aerodynamics 


The problems in the experimentation of rarefied gas are discussed. 


First the extremely large viscous effects in a wind-tunnel nozzle is 


shown. Then the difficulties of flow measurement are surveyed, pointing 


out partioularly the unconventional behavior of the Pitot tube in 
rarefied gas. The performance of a hot-wire anemometer is then studied 
in some detail to show its feasibility. Finally the rules for achieving 


complete flow similarity of rarefied gas flow are formulated. 
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THE BUCKLING 888 SHELLS UNDER 


Th. von Karman and Hsue-shen Tsien 

In two previous papers (Ref. 1 and Ref. 2), the authors 
have discussed in detail the inadequacy of the classical 
theory of thin shells in explaining the buckling phenomenon 
of ere cylindrical alma and tre spherical shells. 
It was show that not only the calculated buckling load is 
4 to 5 times higher than that experimentally observed,but 
the buckling wave pattern found is also different from that 
predicted. it was pointed out, furthermors; that the 
different explainátions for this discrepancy advanced by 
L.H. tonne) ig "ud CM untenable when certain 
conclusions draw from these explainations are compared with 
the experimental facts. The authors are then led ከ7 
ጩጨ the theoretical investigation on spherical shells (Ref.1) 

a model experiment (Ref.2) én thin columns: $y 


with non-linear elastic support/to the belief that the 


by means of tà; non-linear large deflection theory. 
tho tun add 


lates we 
the pepid_Aeepm load necessary 
to keep the shell in equilibrium Once e she started to 


buckle. This characteristic of > ከ load sale st 
of all, that there is a release of elastic M once the 


buckling has started, and thus explains the observed rapidity 


2 


of the buckling process. Furthermore, this characteristic 


also brings in the possibility of a decrease in, buckling 
load when there 5T? slight imperfections in the test specimen 


MÁ 
or when there are vibrations during the testing process. 


In this paper, the authors will show by means of 
an approximate calculation the dropping load characteristic 
of a uniform thin cylindrical shell under axial 


/ 
compression. Consequently, they hope that they have thus 


offered an aceptable explaination of the observed facts. 


Stresses in the Median Surface and the 


Expression for the Total Energy of the 
System 
(4.4) 
If u,v, and ወሃ ድፍ the displacement of & point on the 
median surface of the shell in the axial, ze direction, 
the circumferential, Ys” direction, and the radial direction, 
then the unit strains in the x and y directions, Ex, £y 
and the unit shear ly at a point in the median surface 
can be expressed in the foltowing form, including terms 


up to second order: 
2 


p 775 L ታይ gë 


where R is the radius of the i median mi A of the shal. 
The strdses tenias in the median NC of the shell 


are, however, related to each other by th regsions . 
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^ ማዊ 


ui ውው Sa S A A 


(2) 


= Ae hat ሓመ” Y à Payne nabo. 
a RN oy eA bs cmt been fas 8(1) into Equs.(2) the following 


connections between the crac. the median surface 
PT 


and the wa cud of the median surface are obtained: 


ie sie} Be ዓራ ቀጠ #/ / 
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ibrium conditions 
PET 
of the stresses in the median surfare can be 


expressed by the following equations in case of thin shells: 


RE ቫም” 

Es ds (4) 
72 Ze 

This pair of equations can be guk&emeitesiiy satisfied by 


introducing the well-known Airy stress function, F, defined 
by the relations 


L 2 
E / m 
= ቻ ሃማ Ce ET 
2 / 
Eliminating the variables u and v between Fafis. (3) and (5) 


a relation between the Airy stress function,F and the radial ወሃጮ#። 
medina 


8 ue 
acement “is obtained; 


22 A Ais: DZ Jr - el (945) H TD በረ 


H. 


This equation is very similar to the large deflection 
equation for flat plate derived by the senior author. Ü # €) 
L.H.Donnell (Ref; 5) first obtained this equation in its 

d a Ke 


present form. With a given form of the radial 
Eqw. (6) gives the induced stresses in the medjan surface of 


"pum > extensional 
^ I 


a b 
le-. ET 4 "m 
ክኔ ፈፈ. ጨፊ r ራሪ e A mod ፌሪ 8h 
EC Ner aude ድው calculate ም E elastic energy, it is 


E necessary first to find the expressions for the change of 


nad of curvatures an t twist of the median surface. 

In this paper, the following simplified ae ማች 
| P yee [E E £ 
x ^X, T ow 4 y (4) 
) In Ea. (8), certain additional terms in and F involving 

v are neglected. It was shown by L.H.Donnell (Ref.6) 

d that . hexersfeciar differ from the terms retained 
| in Eq.(8) by a factor He where n is the number of waves in 
i= m- the circumferential direction. For thin cylindrical shells, 
! the value of n is y around 10; theralt the neglection 
is jusitfied. With these expressions for the change of 
the ak T wist d the medien surface, the bending 


Wh for nd መቃ . gesl ጋ 


energ tan be obtained as 


K< ረኛ ET ) 26-0 AL N 


W, 
The lowering of the potential, of the applied force 


on the end of the cylindrical shell can be calculated as 


the produet of the change in length of the shell and the 
lied force. p È 

a im T 
W, =- TO (10) 


The equilibrium condition of the shell is obtained 


by equating the first variation of the sum of énépgies, Di, Wa 


D and Was to zero. In the following calculation, the Rayleigh- T 

Y Ritz method or the energy method ቼ838 be used by assuming a -— 
"i ^ 
|! plausible form for the function w. * 
LE ~ 

4 Calculation of the Total Energy Je 


To obtain a plausible form for w , one has to resort to 


Dr: 
1 the Tesults Seemlexperimanted It is observed that, for large E 
L i 


values of the wave amplitude, the waves are so-called diamond E 


d shaped. 'TRis particular wave shape can be approximately ius 
5 expresséd by 4 SES m 
kI,,ů U.. (1) 

j AR JR 


where the squares are introduced to account for the fact that 


the shell has a definite preference to buckle inward. Eq. / I) 


ይ 


can be rewritten as 
ya 2» 
f Dr ሪሪ #መ M sd emt Fe , (12) 


On the other hand, the classical theory which 48 correct for 
&nfinitesimal vs 1968 of the wave amplitude, requires the 


waves to be of the form 


at ng ni 
= Geo FUR 
| Ld. R (13) WÄ 
* LU 
i Í L. L * : 
| I In order to satiafy this requirement, the wave form assumed የ, |. 
| ` n in the following calculations is of the | 
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where fo, f4, and fo are unknowns to be determined by 
| [ conditions of minimum total Mii a of the system. 
1 fo is introduced inorder to altow the shell to expand 
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The solution of this differential equation can be easily 

ጌ 
btai ይ 2 
obtained as «£i 
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e de 
Fo- ET (TA ME i * ፌዝ om ex] 
D to a, 
ME gt AN 
Using Eg. (5), the stress components in the median surface 


AI de ዖኪ cu 


e ep | ፥”ኞ* pine V R 


PD a , +f 


mY en T E oat ak, 


= G- gels wart , መ= 7732 n GA 


ም ME ቴመ EL r 


C ያ 
4 ET lta ca a f an 


Ty = 101 | 


In all experimental results, the data is usually 


expressed intterm of the average compression stress jin the 


axial direction. It can be easily seen from Eq.(18) that 


#- -፡ 


(39) 


57 ng Using Eq. (3), a Lor the following 


expressions for Ze ana 97 can be obtained: 
DeL H 27 4 
A ፖ/ኳ TIER 


EVN ve Achte. 


By substituting Eds. (14) and (18) into Eq. (20), it is found 
that 


d a NE lA] cosine Terms 


(21) 
#ሃጅ 2 "/ሯ/ HE IT teh) + Coline Term; 


Since y-direesten ¡the circumferential direction, v 
must be a periodic function; therefore the constant te m in 
Ce must beYZero. ጨጨ by 


s0 
Beg [O AGA 


This im the condition s to deterninel. foc ^ 
7 € 
Using Eqs. (18) and (22), the externsional energy of 
the shell is obtained as 
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Using Eqs. (9) and (1&),,the bending energy of the shell 


can be calculated as 


Bm uw | 
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The lowering in potential of the aplied force can be 
obtained by means of Eqs.(10), (18),e#2@1$2and (22): isThe 


result is 


AT. È Sg 2 X 
Zu = - EE ፉ#ሃ፲72 ሥሠ (55) 
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Relation between the Compression Stress 
&nd the Amplitude of Waves 


To find the relation between the compression stress 
and the amplitude ፅያ the waves, the conditions which will 
make the sum of the energies Wy Wo , and Wz a minimum 
have to be obtained. It was found that the calculations 


can be simplified to a certain extend by first est rd 1 . 


SA MEME Aj) T (26) 


This condition will determine toe relation between fo 
m ሪራ de a: 
and fy and fo J 
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Subdituting Eq.(27) into the expressions for W ,W 9 and Ws 
as given by Eas. (25), (24) and (25) and using Ea. (16), 

the total energy d the system 48 expressed finally in 

the following form: 
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The equilibrium conditions are then obtained by differentiating 
the expression for tbtal energy, Eq.(28), with repppct to 
f, and fp , and then set these expressions equal to zero. 
The results can be written in a simpler form by introducing 
the following parameters 


ha 
dll. 


ei dus ecuilibrium conditions are 


palo) | vase t her t 


PR ETE 
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Eliminsting £5 ffom the above eguations, the following 
equation for È is obtained 


Arie A Apt Ao 


where the coefficients are 


TE AN 


(32) 


Thus with a given value of ሥ and # , the coefficients 
for various values of the wave amplitude E can be ftrst 
calculated by using the Eq.(32). Then Eq.(31) can be solved 
for different values of P corresponding to this particular 
value of ana at varbous wave amplitude - When 
the value of ያ is known, Eq.(30) can be used to calculate 


the corresponding value of the compression stress È . 


It is found, however, that the following expression for the 
compression stress which is obtained from Eq. (50) by 


15 
‘eliminating the third powers of 3 ,38 more suitable for 


numerical comput®tion: 


These calculation was carried out for two values of 
the parameter ሥ , the aspect petico of the —waves,—or mone 
defirmittelsy the ratio of the wave lengthscin circumferential 
direction and in axial direction. These values of / are 
1 and 0.5. The former value was chosen, because the 
experiments indicates that at large wave amplitudes, the diamond 
waves havévécual sides. The latter value of ሥ was chosen, 
because the experiments made by N. Nojima and S.Kanemitsu 
(Ref.2) showthe occurance of this type of waves. The 
results of these computations hre shown in Fig.2 and Fig.3 


where the compression stress is plotted against the wave amplitudes. 


Rebel ju + 


> 14 
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The parameter in the figures is 


, 4 . The values written dk 
in the parenthesis is the actual number of waves in circum- ` pr! 


ferential direction for cina 1 = 1000. The dropping 


pi load characteristic is PERE Y evident . JY . if 
A m» det fo top valu ባን Ape A hell de Fut, “ረ Wa 
Zén l i 8 . È Ce a 
Ma ir merde d ar 8 re. ade ab wese er de rd Au 3 VO 
፻ with a cala d YA f^" e Relation between the Compression Stress pi Cé 
Za and the Shortening in the Axial Direction Mo ሥሥ ጨ Ay 
H CM de duo tù gd Mn 
% Although the dropping load characteristic of the 


cylindrical shell is shown very clearly by Figs. 2 and 3, 
the actual load taken by the specimen under a testing machine 
can not be predicted from these figures. In a testing 
machine, the only factor under the control of the operator 
is the distance between the end plates and this is the feestor 
whieksdetarmine geometrical restraining the specimen 

must conf@fm. Therefore, it is clear that the compression 
stress on the shell 1s not determined by the wave amplitude, 
but by the end shortening. The unit end divis icd" 1.0., 


from Eq.(21). It is found that 


| f= S£» IE ibis, 


(34) 


. the aff ci mate labeler put "eS nue ae 
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derd ar ¿noble b do 


Q values of Se and 


4 and 5, the compression stress is plotted against the 

unit end shortenings, for tae P / and 

It is immediately clear from this two figures that if 

the buckling process follows the curves bor once the 

shell starts to buckle the end shorteni 5555 
the tti 


the end NC EE adi This—is fuere 
A = መ ዛ ፉሬ ki ሠረ መ Zar 7 the 


Therefone,-/the shell must first jump to the point P (Fig4 or 
5) which has the same end shortening as the starting point of 


buckling process, but has a much lower compression stress . 


This jump definitely involves a release ስደ elastic energy m 
| | and thus explains the rapidity of buckling process observed Pr 
E | AE ee tig accompanied vibration zho Kees afe ሥራ. | ER 
+ | However, it 48 still difficult to safxthé shell will i A 


A jump Mo a state with narrow waves, as shown by Fig.5, am h m 
1 D 

de [ 

È to a state with almost square waves, as shown by Fig.4.) "m 
4-8 i 


d 

4 Superficially, one might conclude that narrow waves are E D 
gd Yd In. 
174 probabl ecause it gives a muchicompression stress, as M" mY 


can be seen by comparing Fig.4 with Fig.5. But at a 


given end shortening, the true criterion 


EE 


fundamental _concept of mechanics/for the most probable 


SS 


equilibrium position is that the elastic energy stored in 


p 
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shell should be the lowest. An approximute calculation 


of the elastic energy of the shells shows that for values of 

È < 0,012/ s the occurance of narrow waves is improbable, 
because the elastic energy is higher than tnat መ square waves 
at same value of end shortening. However, for F near 
0.6, the elastic energy stored in the shell for narrow waves 
is comparable to that for the square waves at same values of 
the end shortening. B ex indicates sone bap 
narrow waves w m 
fhe—Paiuu—utumweme-putrt—Uf-weibeeg. This is in agreement 
with the experimental findings of N.Nojima and S.Kanemitsu 

as reported in a previous paper by the present authors (Ref.2). 

In any case it is d e that there are equilibrium 

positions of a bucklea cylindrical shell which involve 
much lower average compression eas es that ptocicted 

the beginning of buckling. Take for instance, ei 

E ስ lowest compression stress is given by 


f- 0.194 (È) (85) 


This value corresponds closely to ali pi the experimental 
results obtained by bed E,E.Lundquist and L.H. Donnell. (% ) 
In the two previous papers (Ref. 1 und Ref. 2) the authors 

have pointed out that due to imperfections present in 

test specimen and vibrations in the surroundings, =5= 


preeblo-+ the test buckling load will probably nesrer to 


17 
the minimum load necessary to keep the shell in buckled 


shape rather than the initail buckli load given by the 
classical theory. H x A t ion 
tdi ‘ Ans y ¡mado 


The Effect of the Elastic Characte istic of the 
Testing Machine on the Buckling Behavior 


In the previous paragraphs, it is stated that the 


compression on the specimen is determined by the distance 


between the end plates which is in turn determined by 
| the loading mechanism. However, in actual cases, there is i 
always a certain amount of elasticity in the testing machine 
E and thus the distance between the end plates is not only Pete. 
L determined by the loading mechanism but also by the compression H 
| force wren on the specimen. With the loading mechanism P 
Lee, at a certain position, a compression force on the specimen * 


` PN anfani > 
ሾ Net the specimen; ssuming that the testing machine has 


E i will force the end plates apart and thus reducs tS Gnd shortening 1 È 
a linear elastic behavior, the &ompression load is related FA 
to the end shortening by parallel straight 1ines\for Constant 
» values of, g&y,—number of turns of the loading crank of the T ‘ 
testing machine. If the lading crank is held at a fixed " ጻ, 
position, then the compression load of the specimen must d 


follow the straight line corresponding to this crank position. i 


In Fig.6, the load-end shortening curves of different 
character are superimposed on a system of straight lines representing 
fhe characteristics of thetesting machine. It is evident that 
after the maximum, or initial buckling load is reached, the 
shell will jump to ድክ ila bps state involving much lower 


compression load. The new sa state is determined 


Y Ai 


by two factors: 5 e shape of the load-end shortening 


rated 
the curve B (Fig.6). e elastic character of 


curve. The curve A me will give higher load than 
the testing machine. A more elastic machine will give 
a set of characteristic straight lines having smaller slopes. 


Therefore, in case of Curve A, a more elastic machine will 


result in a higher load, while in case of curve B, a more 1 
elastic machine will result in a lower load. Therefore, 

the buckling behavior of a cylindrical shell 1 determined 

by the elastic characteristics of the shell itself, but also ler. 


by the elastic characteristics of the testing machine. : 8 
se — ob sull 

This fact UC first ripe by the senior author As 

de L. partic 

E.E rr b 


TE-thepe”is 1 5 on n en, 


Dn mim 


Concludi Remarks 


In the previous paragraphs, the authors have shown 
that there are equilibrium positions in the buckled shape 
invólving much lower load than the buckling load predicted by 
classical thoeyy, and thus if the specimen is slightly 


imperfect, it is reasonable to expect much lower buckle loade 


ጨፉሖ a | a i 
“ef Ythe order served in the laboratory. They have 


also pointed out that the elastic behavior of the testing 
machine has a profound influence on the buckling phenomenon 
and this might be another cause of the large scattering of 

e L . obtained by different experimenters. 
However, due to complexity of the problem, the results 

given in this paper can be only considered as a rough 
approximation and most of the F ped a made 46 qualitative 
rather than quantitative. To put the new theory on a 
solid footing, a more accurate solution of the differential 
ecuations of equilibrium is necessary. Particular attention 


must be given to the calculation of the elastic energy stored 


“ጋሪ in the shell, because it is found that £imel*factor which 


4 ኣ determines the modt provable equilibrium DEE be 


| the compression stress on the specimen nor the wave amplitude of 
maths. / 


the buckled cylindrical shell, but rather the elastic energy 


stored. a ጋና "YI a fe 


Perkeps, a meme inquisitive tmrestigetor E e pleased 
by & rigorous proof of the validity of all the large deflection 
equations. equations weed are established by 
ei intuitive -arguments rather by systematic reasoning. 

For instance, it is not certain whether the curvature of the 
shell has to be calculated more accurately by taking into account 
the second order terms, or the extensions of median surface 
should be more accurately determined. It is the belief of 
the puesert authors that an investigation of these problems 
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by starting from the general non-linear elasticity theory 
kag Gosserat (Ref.7) and others vide very EA 

The senior author has already expressed this part 

opinion in his 1949 Gibhs Lecture Age: 1 % dean. sla auf 
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THE BUCKLING OF THIN CYLINDRICAL SHELLS 


UNDER AXIAL COMPRESSION 


Theodore von Karman and Hsue-shen Tsien 
California Institute of Technology 


In two previous papers (Ref. 1 and Ref. 2), the authors have discussed 
= detail the inadequaoy of the classical theory of thin shells in explain- 
ing the buckling phenomenon of cylindrical and spherical shells. It was 
Shown that not only the calculated buckling load is 3 to 5 times higher 
than that found by experiments, but the observed wave pattern of the 
buckled shell is also different from that predicted. Furthermore, it was 
pointed out that the different explanations for this discrepancy advanced 
by L. H. Donnell (Ref. 3) and W. Flügge (Ref. 4) are untenable when certain 
conclusions drawn from these explanations are compared with the experimental 
facts. he- uuthore-are- Chen led By & theoretical investigation on spherical 

d Us autioro ure iek na 
shells (Ref. 1)/to the peated that e buckling phenomenon of curved shells 
can only be explained by means fot non-linear large deflection theory. This 


N 


AA! alenke 
point of view is—further substantiated by 4 model experimentson bn columns d 
TTT dee EE E 
with non-linear elastic support (Ref. 2). Itis-evident—frem_these—two 
invostigations—bhet the load necessary to keep the shell in equilibrium 
drops very rapidly with increase in wave amplitude once the shell started 
to buckle. This characteristio shows, first of ell, that there is a release 
of the elastic energy stored in the shell once the buckling has started, 
and thus explains the observed rapidity of the buokling process. Furthermore, 
it also brings in the possibility of a decrease in the buckling load when 
there are slight imperfections in the test specimen and when there are 
vibrations during the testing process. 


In this paper, the authors will show by means of an approximate 


= 


calculation that in case of a thin uniform cylindrical shell under axial 


compression, the load sustained by the shell drops with increasing deflection. 
Consequently, they hope that they have thus offered an acceptable explanation 


of the observed facts. 


Stresses in the Median Surface and the 
Expression for the Total Energy of the System 
Let X and / be measured in the axial and the circumferential direction 


in the median surface of the undeformed cylindrical shell and 4 , Y , and 


W be the components of displacement of a point on the median surface of e 


the shell in the XV direction, the / -direction, and the radial direction. b 


Then the unit strains in the X end y reete 2 , > and the mit e 


shear A &t & point in the median surface can be expressed in the following 


forms, including Mer D 2 ደው order: 


JE 289 0) 
ፉ QU 255 
hey = ዳም * 2x è; 
where R is the AM of vis deformed median surface of the shell. The 
stresses and the strains in the median surface of the shell are, however, 


related to each other by the following equations: 


(2) 


= 222 


where Æ is Young's modulus of elasticity and Y is Poisson's ratio. There- 


fore by substituting Eq. (1) into Eq. (2), the following connections between ' 


the components of stress in the median surface and the components of dis- 


placement of the median surface are obtained: I 


=> 

AE 
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It is generally accepted that the equilibrium conditions of the 


(3) 


Stresses in the median surface of a thin shell can be approximately 
expressed by the following equations: 


26 , 22. 
or T ia 


2 07 
ER + daa 
This pair of equations can be satisfied by introducing the well-lknowm 


(4) 


Airy's stress function, Fk, 72 defined by the relations 


LE "E y 
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Eliminating the variables ፊ and // in Eqs. (3) and (5) a relation between 
Airy's stress function Fx, 7 and the radial component of the displace- 
ment, j/ is obtained: 

PERE „„ EE KEY 

BIZ, wy ፖ Ge # KR ia RIL (8) 
This equation expresses the condition of compatibility between stress and 
strain. When P, it reduces to the corresponding equation for a flat 
plate derived by the senior author (Ref. 5). L. H. Donnell (Ref. 5) first 
obtained this equation in its present form. With a given form of the 
radial component of the displacement, i/ , Eq. (6) gives the induced stresses 
in the median surface of the shell. 

For one complete wave panel, the extensional elastic energy ZA corre- 


sponding to these stresses can be written as 


a,b 
Me GE lte: ue . e (7) 
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where and / are the half wave lengths in the axial and the oiroumfer- 


ential directions respectively. 

To calculate the bending elastic energy, it is necessary to find the 
expressions for the change of curvatures and the unit twist of the median 
surface, In this paper, the following simplified ea will be used: 

VE EZ 
KE dy i, Ch A4 (8) 
In Eq. (8), certain additional La in a" involving V are 


Y, 
neglected. It was shown by L. H. Donnell (Ref. 6) that those negleoted 


terms differ from the terms retained in Eq. (8) by a factor # » where 
J is the number of waves in the circumferential direction. For thin 
cylindrical shells, the value of H is around 10; therefore the neglection 
is justified. With these expressions for the change of curvatures and 
the mit twist of the median surface, the bending energy # for one 


complete wave panel ይ ye written as 


RE RE ያ 


The lowering of the potential # of the applied force on the end 
of the cylindrical shell can be calculated as the product of the change in 
length of the shell and the applied force. Therefore the following 


expression is obtained for one complete wave past 


M, = -44 f) 4 [324 (10) 


The equilibrium condition of the shell can be obtained either by 
equating the first variation of the sum of the energies ZA » ZA and # to 
zero, or by actually analyzing the moments and the stresses in the median 
surface of the shell. Using the approximations stated previously, Donnell 


(Ref. 6) derived the equilibrium equation &s 
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where ፆ is the external radial pressure on the surface of the shell. In 


the oase concerned, pro » then using Eq. (5), the second equation 


connecting Airy's stress function Fay ያር the radial component of dis- 
placement, A is obtained as 


EEE fije bz dr" ETES “ሪያ ያ ያ. 


When P, Eq. (12) reduces to the corresponding equation for a flat 
plate. 


There are two different ways to solve the problem of the buckling of 


a thin uniform cylindrical shell under axial compression. The more exact 
| method is to solve Eqs. (6) and (12) simultaneously, using appropriate l * 
d boundary conditions. The approximate method is to first assume a plausible ! ei i o 
| function for , with undetermined parameters and then use Eq. (6) to |1፪ 
determine the stresses in the median surface of the shell. The energies pm 
E ! MW» W, end JY; can then be 6819918594 ty means of Eqs. (7), (9), and GITA 
c (10). The undetermined parameters can be ascertained by the condition that 
| the sum of energies N, » A , and LA must be a minimum. This approximate 


} method will be used in the following calculations. ^ 


gi Calculation of the Total Energy 


A d 
8 To obtain a plausible form for W , one has to resort to the experi- ) 
] D L mental results. It is observed that, for large values of the wave amplitude, Mk 
ul | i 
Í the waves show a so-called diamond shaped pattern. This particular wave ] TN 


shape can be approximately expressed by å 


^ Mon ep IN Um p. 
dh * Cos XU መ (35) 


Ko 


=6= 
where the squares are introduced to &ccount for the fact that the shell 


has a definite preference to buckle inward. Eq. (13) can be re-written as 
t NC Wu. TUE S, A E go A 
R IE nM + pende us E] 


On the other hand, the classical theory which is correct for infinitesimal 


(14) 


values of the wave amplitude, requires the wave to be of the form 

A GEET 

R = (49 id la የ (15) 
In order to satisfy this requirement, the wave form assumed in the following 
calculation is 


# 
x. (A+ 05 Zen tg hoz) lot iz) 
where # " Á , 4 are unknowns to be determined by conditions of 
minimum total energy of the system. Á is introduced in order to allow 
the shell to expand radially. The amplitude of the wave pattern defined 
as the maximum difference in the radial deflection, / is evidently given 
by £ . The wave lengths in the axial and the circumferential direction 
are ፊደ and መ respectively. Hence the number of waves &long the 
circumference of the shell is equal to 7, . It is evident that no end 
effect can be &ccounted for by this form of wave pattern, and therefore 
the following calculation really corresponds to the case of a very long 
cylindrical shell. This simplification is justified by the experimental 
findings of N. Nojima and S. Kanemitsu as reported in a previous paper 
(Ref. 2). It was found that there is no appreciable length effect when 
the length of the cylindrical is greater than 1.5 times the radius of the 
shell. Furthermore, it is seen that by setting =/=, Eq. (16) is 


reduced to Eq. (14); while by setting Z. Hes and 2 Ti =0 » Eq. 


(16) is reduced to Eq. (15). With other values of these parameters, wave 
patterns intermediate between these two limits can be obtained. 

Substituting Eq. (16) into Eq. (6), the differential equation for 
Airy's stress funotion Fs, y) is obtained: 


arfa) =- Pr EED Cor ca Zt 
+0 cnp T Gn Won DE +H caz TEN daz 


where A. 4 , the aspect ratio of the waves. If , the waves are 
longer in the circumferential direction; if ሥዶሩ/ » the waves are longer 


in the axial direction. The coefficients in Eq. (17) are given by the 


following relations: 
ል- pAn- (3H 4) 
2 2 
B= j4 » 
o- 34m lH A) - 14 
v- ፋሪ (4444) 
o= LA ry 
and E= n° (24 +h) 
The solution of this differential equation can be easily obtained as 
ምሃ A ረ Po, c. 7 
e EP TE en e aL + rg መሻ 


2 mr, N 221% % IN coe PM 19) 
"ቃፖ > nH + 27 E ሥሯ "iech die. 
EPA 2 
2 2 f. 
Using Eq. (5), the stress components in the medien surface can be written 


85 > D 

OE ken ku. EE fe E 
#ሪ 

en Men EET 


(50) 


OCH a Reet oye መም cnt 
+ 2 ee za 


y- Ko or . 
E. MEM EF LL tim e cia 2E] 


In all experimental work, the data is usually expressed in terms of 


the average compression stress DY in the axial direction. It can be easily 
seen from Eq. (20) that 


i (21) 


Using Eq. (3), the following expressions for 24 and 2 can be 
obtained: 7 


ee, 
Lb bu “ማም iù 


By substituting Eqs. (16) and (20) into Eq. (22), it is found that 


4ፉ = ve 7321/11 ^ gc AA * Terms of periodio 


funotions 


(28) 
XX giat 4 E). pur n ከ ር በል 


Sinoe is measured along the circumference of the shell, A must be a 
periodic function of / ; therefore, the constent term in # must be 


equal to zero. ዐዮ 7 


cata PAGA +744 l 0 (24) 
This is the condition for the determination of X 


Using Eqs. (7), (20) and (24), the extensional energy # of the shell 


is obtained as 


pu — 4/6 ME ም dl) de 
-Inl GA + ALA HE] . 4 SS (25) 


^g t “መያ” o] 


Using Eqs. (9) and (16), the bending energy LA of the shell can be 
calculated as 


ga cir et o 


The Lewering-of_ihe potential of the applied force oan be obtained by 
means of Eqs. (10), (20), (23) and (24). The result is 


iud በፈሪ ርዕ E ap + P a „ 
nw EF). Wi 


Relation between the Compression Stress and 
the Amplitude of Weves 
To find the relation between the compression stress and the amplitude 
of the waves, the conditions which will make the sum of the energies VAI " 
W, » ena jf, a minimum have to be obtained. It was found that the 
calculations can be simplified to & certain extent by first using the 
condition that the sum of energies must be minimum with respect to Z . 


EMM M;)=0 


This condition determines a relation between Pd and 4 and A » which 


Or 


(28) 


can be written &s: 


f, Gt E44 + 


Using this relation and Eq. (24), it is easily seen that =. In other 
words, the shell will expand radially to such an extent that the average 
of the circumferential stress fe is equal to zero. Substituting Eq. (29) 
into the expressions for # " # and LA as given by Eqs. (25), (26) and 
(27) and using Eq. (18), the total energy of the system is expressed finally 


Se doct 
SC MEV መገሻ A A fer 
+ #ሯራሃ zeit RE dën ER 
thes tt ሄቻ * È Ley] “#ሯ ይታ ICT] -— 
wf lh d ge laa EL 
ETE ^ 


«a o PR ያ + ttg] 


The equilibrium conditions are then obtained by differentiating the > ላ 
expression for total energy, Eq. (30), with respect to X and 4 , and 


then set these expressions equal to zero. The results oan be written in 


TNT à a simpler form by introducing the following parameters: 
8): ,) 
i * 
| A i$ ያ 
í መ  ኢ= = A = EM A 
4 „ T A E (31) 


where $ is the wave amplitude of the buckled shape of the cylindrical 


shell. Then the equilibrium conditions are 


Helle i- Bez 2 /#ይ- Go Sale 

1% ታያ CH Zeie, ## ጋሪ 2 NI 
m Te 2 Lie +t tex ff EDIT f ff 

* ae V [0g + Habeo Ma M (2) 
By ‘Up 4) = GET eg 2 ES +e 


መሪ 2 Tianhe 
AE “tigen ll #/#፡ ታች 


ltd xa JU (pos d ont] 


Eliminating E from the above equations, the following equation for TE 
is obtained 


ASTAS e 


p E : « 
NIE ën el 
ER - ME fr) 


n (UI ug GE 2 “2 ZA Bk, ረጋ Y: tane 
A, OE Fp does Tif Fi Lig # 
“#፡ጋ 7 zi s E oof 


Therefore, when = 


(55) 


„ 1. e., when the wave amplitude approaches zero 
Eq. (32) gives 4,=/-0 ; ena 
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Substituting into Eq. (51) it ie seen that SÍ. or £- #4 P 
Putting this relation between £ and 4 into Eq. (14), the wave pattern is 
reduced to that represented by Eq. (15), i.e., the classical wave pattern 
for infinitesimal wave amplitude. 

With a given value of A end 7 » the coefficients for various values 
of the wave amplitude Z can be first calculated by using the Eq. (34). 
Then Eq. (33) can be solved for ም corresponding to this particular set of 
values of ሥ and 7 at various wave amplitude f . When the value of 
is known, Eq. (32) can be used to calculate the corresponding value of the 
compression stress £g . It is found, however, that the following 
expression for the compression stress which is obtained from Eq. (30) 
by eliminating the third powers of $e is more suitable for numerical 


computations: 


RI a a 

TER digo dr fe 
re ln us) mi e n 
MA iz tp) 


Mil- 


3 
Therefore, when ያቃላሪ , i.e., when the wave amplitude becomes very small, 
Eq. (35) reduces to 


f Ee 2773827 Gel Se (56) 


The minimum value of 2 average compression stress is given by 


NE 
del Wer (27) 


which is the well-known Si from the classical theory of infinitesimal 
deflections. This minimum value is obtained when 


2)2 

p = ሪቱ Hi-v) (38) 
It is interesting to notice that for infinitesimal wave amplitude, the 
minimum value of average compression stress is not determined by separate 
parameters ያመ , but by a combined parameter shown in Eq. (38). 

These calculations were carried out for two values of the parameter 
pe » the ratio of the wave lengths in circumferential direction and in 
axial direction. These values of A are and 05 . The value / 
was chosen because the experiments indicate that at large values of wave 
amplitude, the diamond waves have almost equal sides. The value 0% 
was chosen to investigate the possibility of occurrence of narrow waves. 


The results of these computations are shom in Fig. 2 and Fig. 3, where 


the compression stress qb is plotted against the wave amplitudes 2 . 


The parameter in the figures is » . The values written in the parenthesis 


after 7 is the actual number of waves % in ciroumferential direction for 

R/£ = 1000. For a given value of odp, i.e., a fixed size of the wave, 

the load sustained by the shell, ያ first decreases as the wave amplitude, 
E , is increased, After a minimum is reached, the load will rise with 


increase in weve amplitude. When the waves are larger, the initial buckling 


Me iie 
ý d load, i.e., the value of # at 7=0 , is higher. However, the minimum 
E? load reached tends to a lower value, except for /< 0/67 and AO. For 
7 E 4 (“0.5 » the lowest value of the minimum load is not yet reached at ge ec, 
E 
{ da i The Relation between the Compression Stress 


and the Shortening in the Axial Direction 


Although the load characteristic of the cylindrical shell shown in 

the Figs. 2 and 3 give the possible equilibrium relations between load and 
amplitude of the deflection wave, the actual behaviour of a specimen ina 
testing machine cannot be directly seen from these figures. Ina testing 
machine, the only factor under the control of the operator is the distance 
between the end plates; this is the geometrical restraint the specimen must 
conform. Therefore, in order to determine the behaviour of the specimen the 
compression stress will have to be plotted as function of the end shortening. 
The unit end shortening, £ , i.e., the total shortening in one wave length 
of the shell in axial direction divided by the wave length, can be easily 


calculated from Eq. (23). It is found that 
ER ድያ 2 2 3 7 
„ R4 (59) 
This equation for the unit end shortening contains only quantities already 
found such as the values of f£ and £- . In Figs. 4 end 5, the compression 
ደሮ x 
stress El is plotted against the unit end shortening "2 ደው / and 
| ሥራሥራሪ „ respectively. It is immediately clear from this two figures that 
if the buckling process follows the curves draw, after the shell starts to 
buckle the end shortening has to decrease. In other words, the end plates 


of the testing machine have to move apart. Therefore, the process of 


buckling in this region is highly unstable, since before the operator has 


—— del A 
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times to separate the end plates, the shell will jump to the point Z 


(Figs. 4 or 5) which has the same end Shortening as the starting point of 


the buckling process, but has a much lower compression stress. This jump 
in equilibrium positions involves a release of elastic energy and thus 
explains the rapidity of the buokling process observed and the accompanied 
vibration. 

However, it is still difficult to s&y whether the shell will jump to 
a state with almost square waves, as shown in Fig. 4, or to a state with 
narrow waves, as shown in Fig. 5. Superficially, one might conclude that 
the narrow waves are the more probable shape because it gives a mich lower 
compression stress, as can be seen by comparing Figs. 4 and 5. But, 


according to the fundamental concept of mechanics, the true oriterion for 


the most probable equilibrium position compatible with a given value of the 
end shortening is that the elastic energy stored in the shell should be the 
lowest. An approximate calculation of the elastic energy of the shells 
shows that for values of Tea ll , the occurrence of the narrow waves is 
improbable, because the elastic energy stored is higher than that for the 
Square waves &t same value of the end shortening. However, for ያ near 
0,6, and VEIL, / the elastic energy stored in the shell for the narrow waves 
is comparable to that for the square waves at the same value of the end 
shortening. This indicates the possibility of the appearance of narrow 
waves during the very initial stages of the buckling process. 

In any case, it is certain that there are equilibrium positions of a 
buckled cylindrical shell which involve much lower average compression 
stress B than that &t the beginning of buckling. For instance, in the 


case of square wave pattern, the lowest compression stress is given by 


L 0/94 Z 


E (40) 


18 
Incidentally, this value corresponds olosely to most of the experimental 


results obtained by L. H. Donnell (Ref. 3) and E. E. Lundquist (Ref. 7). 


The corresponding value of the parameter 7 which determines the 


number of waves is equal to 0.225. In case of 22 = 1000, the number 
of waves, JL , will be 15 which also agrees well with the experimental 
evidence, For this particular value of the radius to thickness ratio, the 
number of waves along the circumference decreases from the JM = 26 at the 
beginning of the buckling to 7 = 15 at the calculated minimum buckling 
stress, This gradual increase in the size of waves with the unit shortening 
is also observed by the experiments reported in an earlier paper (Ref. 2). 
It is particularly interesting, however, to trace the gradual change 


in the wave pattern during the buckling process. Figs. 6 and 7 show the 


“መ lines of equal deflection of the wave surfaces corresponding to different 
equilibrium states for two values of the aspect ratio of the wave pattern, I 
pelo and fees - These particular equilibrium states are denoted in 
Figs. 2, 3, 4 and 5 by a small circle in order to indicate their relative 
position during the buckling process. It is seen that there is a rapid 
shift from the rectangular waves bounded by lines, ፡ = const., and E = 
const., as predicted by the classical theory for infinitesimal wave ampli- 
tudes, to staggered rows of circular or elliptical waves. Whereas, the 
rectangular waves are directed alternatively inward and outward, the circular 
or elliptical waves are all directed inwards. The transition is practically 
M completed for = 0r £, i.e., when the wave amplitude is only 4 or 6 times 
ያ the thickness of the shell. The occurrence of such inwardly directed 
[ ከ] circular and elliptical waves at this stage of the buckling process is in | 


good agreement with the experimental observations (Ref. 2). If the experi- 


" ment (Fg) is continued to larger deflections, these staggered waves ዘ 
be on | We A | 
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obtain characteristic diamond shapes, The present approximate theory fails 


to give these sharp diamond shaped waves. It is obviously not sufficiently 
exact for such large deflections, Furthermore, when these diamond shaped 
waves occur, the load on the specimen actually falls to a very low value 
such as ድ = 0.06, whereas the theory shows a slight increase of the 
stress at least for the case ሥራ = 1.0. Therefore, the present calculation 
can be only considered as a good approximation to the earlier stages of 
buckling when the wave amplitude is only a few times the thickness of the 
Shell. Nevertheless, it reproduces the characteristic features of the 


buckling process observed in the laboratory. 


The Effect of the Elastic Characteristio of 


the Testing Machine on the Buckling Phenomenon ` m i 


It was stated in the previous paragraph that the state of the specimen 


is determined by the distance between the end plate and that this distance 


is the independent parameter controlled by the experimenter. This statement TEN 
is correct only insofar the elasticity in the mechanism of the testing d ቁ d 
machine is neglected. Since there is always a certain amount of elastic gn 
deflection in the loading mechanism and this deflection is a function of M 
the load. Hence, if, for example, the loading crank is held at a certain o | 


position, the compression force acting on the specimen will foroe the end 


plates apart and thus reduce the amount of end shortening of the specimen. 
+ | The actual shortening is determined by the load-deflection characteristics 
| of the specimen &nd the testing machine. Assuming thet the testing machine ‘a 
A has a linear elastic characteristic, the compression load is related to the 
| end shortening by parallel straight lines, each iine corresponding to, say, ብ 


a constant number of turns of the loading crank. If the loading crank of d um 


——— ———— 
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the machine is held at a fixed position, corresponding values of the 
compression load and end shortening of the specimen must lie on the straight 
line for this orank position. If the load-end shortening characteristics 

of the specimen itself &re given, it is evident that the equilibrium positions 
of the entire system are determined by intersections of the curves re- 
presenting the characteristics of the specimen with the straight lines 
representing the characteristios of the machine. 

Fig. 8 shows representative curves for the characteristics of the 
Specimen and two families of straight lines representing the characteristics 
of two différent testing machines. It is evident that after the maximum or 
initial buckling load is reached, the shell will jump to a new equilibrium 
position involving much lower compression load, But this new equilibrium 
position is determined not only by the load-end shortening relationship and 
also by the elastic characteristic of the testing machine. A more elastic 
machine will give a set of characteristic straight lines with smaller slopes. 
Therefore, in case of curve A (Fig. 8), & more elastic machine will make 
the shell to jump to & higher load, while in 0886 of curve B, & more elastic 
machine will make the shell to jump to a lower load. This influence of the 
elasticity of the testing machine has been discussed by the senior author 


in connection with the plastic buckling of colums (Ref. 8). 


Conoluding Remarks 
In the previous paragraphs, the authors have shown that there are 


equilibrium positions in the buckled shape involving much lower load than 
the buckling load predicted by the classical theory, and thus if the 
specimen is slightly imperfect, it is reasonable to expect much lower 


buckling loads. They have also pointed out that the elastio characteristic 


=19- 


of the testing machine might have quite a large influence on the buckling 


process and this might be another cause of the large scattering of the 

data obtained by different experimenters, However, due to the complexity 
of the problem, the results given in this paper can be only considered as 

& rough approximation and most of the statements made are qualitative 
rather than quantitative. To put the new theory on & solid footing, & more 
accurate solution of the differential equations of equilibrium is necessary. 
Particular attention must be given to the calculation of the elastic energy 
Stored in the shell, because it is found that the most probable equilibrium 


depends on the magnitude of the elastio energy stored in the various 


equilibrium positions compatible with the constraint exerted by the loading 
process. 
d Furthermore, an inquisitive mind will, perhaps, be pleased by & 


rigorous proof of the validity of all the large deflection equations. 


M LM 

y 

U These equations are established by intuitive arguments, not by systematic CAI 

7 reasoning. For instance, due to the appearance of sharp curvatures in the RE 
Jo 


diamond shaped wave surfaces at large deflections, it is not certain whether 
the curvature of the shell has to be calculated more accurately by taking 

into &coount the second order terms, or the extensions of the median surface 
Should be more accurately determined. It is the belief of the authors that 


an investigation of these problems by starting from the general non-linear 


theory of elasticity developed by G. Kirchhoff, J. Boussiwesq and others 

M is very desirable. The recent work ከሃ R. Kappus (Ref. 9) is a noteworthy 

C contribution in this field of investigation. The senior author has already 

expressed this opinion in his 1939 Gibbs Lecture (Ref. 10) given before the * 


American Mathematical Society. 
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Section 6 


Take-off from Satellite Orbit 
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Take-Off From Satellite Orbit 
H. S. Tsien* 


Danid wd Mene ያሄ...ፈ.... 23 Edi, Mos. LE RL 


ummary 


The mass ratio or the characteristic velocity for the take- off of 
a space ship from the satellite orbit is computed for two cases: the 
radial thrust and the circumferential thrust. The circumferential 
thrust is much more efficient in that the required mass ratio is much 
less than for the radial thrust. Both cases show, however, an increase 
of the required mass ratio and the characteristic velocity with a 
reduction in acceleration. With circumferential thrust, the character- 
istic velocity increases by a factor of two when the acceleration is 
reduced from 1/2 g to 1/3000 g. 
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For take-off of a rocket from the earth surface, it is convenient to 
have the initial trajectory in the vertical direction and then the thrust should 
be considerably larger than the Ínitial weight of the rocket. Depending upon 
the relative magnitudes of the aerodynamic drag and the weight, the initial 
ratio of the thrust and the weight should be between 2 and 3 for minimum 
expenditure of the propellant. The situation is quite different for a space 
ship taking off from the satellite orbit: Ina satellite orbit, the gravitational 
attraction is completely balanced by the centrifugal force, and the vehicle is 
effectively in a weightless state. This fact has led many fanciers of inter- 
planetary travel to conclude that take-off from satellite orbit requires only 
a very minute thrust. For instance, L. Spitzer (Ref. 1) proposed a nuclear 
power plant for bis space ship to be accelerated at only 1/3000 g. Another 


example is the extensive discussion of interorbital transport technigues by 


H. Preston ~ Thomas (Ref. 2) based upon the assumption of egually small 


acceleration. On the other hand, W. von Braun (Ref. 3) seems to prefer a 
very much larger acceleration of approximately 1/2 g for take-off from the 
satellite orbit. 

The magnitude of the acceleration has a strong bearing on the optimum 
type of power plant to be used: The ion-beam rocket is only possible for 
very smallacceleration, while for moderate acceleration, chemical rocket 
is required. Therefore the question of the magnitude of acceleration is an 
important one for interplanetary flight. The purpose of this note is to 
8 the relation between the acceleration and the mass ratio 
required for escape from the earth's gravitational field, starting from the 
satellite orbit. It is hoped that the present investigation will give the future 


generation of astronautical engineers a mere rational basis for designing 


space ships. 


Basic Equations 
The problem considered is the motion of a space ship under the influence 


of the rocket thrust and the gravitational attraction of a single massive body, 
say the earth. Then if the rocket thrust is in the plane of trajectory, the 
trajectory of the space ship will remain in a plane. Let the position of the 
ship at any time instant t be given by the polar coordinates r and f : r 
the distance from the center of attraction, and Ü the angular position. If 
the components of the rocket thrust per unit mass of the vehicle are R in the 
radial direction and © in the circumferential direction, and if g is the 


magnitude of gravitational attraction at the starting satellite orbit 7 = 7; 


(Fig. 1), then the equations of motion of the space ship are 


ለ) 


#/"1/7” 70 a 


By using the subscript o to indicate quantities at the starting instant 
t = 0, the equilibrium condition of the satellite orbit is given by 


CM 


Initially, the radial velocity is zero, i.e., 


0% =p (4) 


These are the initial A: or the space ship to have sufficient energy 
to escape the earth gravitational field at the end of the powered flight, the 


sum of the kinetic energy and potential energy must vanish at the end of the 


-3- 


end of the accelerating period. Let that instant be denoted by the subscript 1. 


Thusat t- ti 
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With any specified variation of the thrust forces R and © as functions 
of time, the above system of equations determine completely the take-off 
trajectory of the space ship. In the following sections, two special cases of 
practical significance will be discussed in detail: the case R - constant, 
© #0, purely radial thrust; and the case R < 0, © = constant, purely 
circumferential thrust. 


Radial Thrust 
If the thrust is always radial and is propertianal to the instantaneous 


mass of the vehicle, a non: dimensions l thrust factor P d can be introduced as 
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Eq. (8) can be immediately integrated and by using the initial condition of 


Eq. (3), the result of integration is 
al - ሠ)” 


ir 


By substituting this equation into Eq. (7), the final equation for ያ 18 
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p al 75 75 

The nondimensional radial velocity is d PI dt. This is related to the 
physical radial velocity dY/dt as follows 
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Eq. (40) can be rewritten as 


Ld ያር POPE 
##፡2።ሥ e e 
Since de /dt - 0 when t - 0 and $= 1 according to Eq. (4), the result of 


integrating the above equation is 


t- Uh) 3-4) E 


Therefore the non-dimensional time * can be calculated as a function of the 
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radius e as follows 
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With Eas. (9) and (55), the end condition of Eq. (5) can be written as 
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Or simply 
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Then the velocities at the end of acceleration period are 
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The time ጣ for the powered flight can be obtained from Eq. (35) by setting 
the upper limit of integration to f}. Thus A nodi / me sok prt SR 


ግ “ም 22 e (a8) 


where F and E are the elliptical integrals of first kind and second kind 
respectively. ® 
If M(t) is the instantaneous mass of the space ship, and c the effective 


exhaust velocity of the rocket, then 


„ 


Therefore the mass ratio M,/ Mi can be calculated ag follows: 
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By using the result of Eq. (16), 
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When the acceleration is very large, per 1, the integrand in Eq. (38) can 


bè expanded in terms of this parameter, Then the mass ratio is calelaled as 


7 


ያ ሣን uL Ce (38) 


በሪ |ፃ 
The relation of Eqs. (ትቸ) and (18) is plotted in Fig. 2. For firs ራ i 


the mass ratio becomes infinite. The reason is that at this value oá 


* The author is indebted to Dr. Y. T. Wu who kindly supplied the relation 


of Eq. n 
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acceleration, there is a radial position where the thrust force is equal to 
the gravitational attraction and one further increase in the energy of the 


Therefore the radial thrust per unit mass, if maintained 
With 


vehicle can occur. 
constant throughout the powered flight, should be larger than 1/8 g. 
increasing thrust, the required mass ratio for escape from the earth's 
gravitational field. decreases. VThe asymptotic value of Join, A ) is pre . 

However there is no appreciable improvement- in going t to higher thrust than 
n 


This 2 dependence of the mass ratio upon the #ር6 616 28636 facto 


Circumferential Thrust 
If the thrust is always circumferential and proportional to the mass of 
ኤቲ, 4“ anek Hat 
L1? 


the vehicle, then + mw Md ሥሠ Y tai 
Ge ቾች (19) 
By using the same non-dimensional variables as defined in Eq. (6), the 
al 
(59) 


equations of motion are 


Te 149-5 
(21) 
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The initial conditions of Eqs. (3) and (4) are 
[77 G= l, TB (22) 
a 
(23) 
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Therefore Eq. (2Q) gives another initial condition that 
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By eliminating £ from Eqs. (29) and (23), 


Lloro 


This is a third ፈር differential equation with three initial conditions 

specified by Eqs. (23) and (24). No simple general solution can, however, 
will de tom tend Kal 

be obtained. The following discussion is-thus centered around approximations 
do 

that are valid for large values of y or amall values of y . 

For very large values of Y , the acceleration period is expected to 
be short and the change of the radial position to be small. Then the value of , 
must be very close to the initial value of unity. By taking e to be unity, 


DCH 


Eq. ፎኔ becomes 
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n 
$where C is the itegration constant. C, however, must be l because of 
the initial condition of Eq. (23). The appropriate approximate solution for 


p for very large y is thus 
b 
M ድ ሰች uut A? 
OS TAL ቃይ ቻቹ . (25) 
To obtain higher terms in this power series, the usual series substitution 
method may be used, The calculation is somewhat lengthy and therefore 
will not be reproduced here. The result is 
sc? 


= iu KE ta mp an (26) 


By using the result of Eq. (26), the radial velocity is obtained by 
differentiation. Then Eq. (28) gives the circumferential velocity. The end 
x condition of Eq. (5) can be modified into the following more convenient form 
by multiplying it by]? 


o- EER ae] 


By substituting the solution of Eg. (28) into this condition, an eguation for 
determining 7, is obtained: 
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The mass ratio M, /M, can be calculated in the same way as in the 
previous section and can be determined through the new parameter x defined as 


follows: 


an 


ips s Ala M)= Vn = ቻ Ger» 


Eq. E then can be written as 
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Since the calculation is designed for large values of V , the appropriate 
ሙሪ 
expansion of x'should be a series in inverse powers of ሃ . Eq. (29) 


suggests specifically 
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where and x are constants independent of V . By 


substituting Eq. (38) into Eq. (39) and equating equal powers of Y , the 


following set of equations results: 
su 
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The explicit numerical solutions are then 
m 
X = #5 ፦/ = ፈቆ/#ኔ/ 
x” = 0001349 
x= -0.000047 


This completes the calculation of mass ratio for large values of the 
acceleration factor V . 


For the other extreme case of very small values of Y , itis to be 


expected that fara. (A) the acceleration will be very small, andthe tArm 


WW RS will be very much smaller than L . Therefore a good approximation 
> 


of Eq. (24) at small y is 
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The solution of this equation with the initial condition of f =lat Z 0 is 
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Therefore 


3 7 
5 | ሥህያ (36) 


At 7-0 , the radial velocity and the radial acceleration are thus not zero 
as required by the initial conditions of Eqs. (23) and ሠኔ They are however 
very small, because Y is very small. Therefore the solution of Eq. (3%) 
is a good approximation to the exact solution. 

To the same approximation, Eq. (20) becomes 
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y X ን (n 


CE - Small — 
This means thatipecause of the extremely siew acceleration, the T 


2 
\ force per unit mass 7 47 practically balances the gravitational 
attraction at every instant) The end condition of Eq. (5) can then be written 


34 
(38) 
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where X is VI, . The appropriate solution for X is then 
A 


X= /፦ (ayy 


Since the mass ratio 22 , is related to X by Eq. (28), Eq. dai 
actually gives the mass ratio for escaping the gravitational field with very 
small acceleration. 

The የ ረ plotted against Y in Fig. 3, using both Eq. (30) 
with Eq. (34) and Eq. ($9). When Y approaches zero, X approaches 1. 
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When ሃ is very large, X approaches WEST T S Y increases, X and 
hence the mass ratio MM, decrease monotonically. Therefore, same as 
the result for purely radial thrust, there is a strong influence of the magnitude 
of acceleration on the required mass ratio. However as far as decreasing 
the mass ratio is concerned, there is no appreciable advantage in using Y 
greater than 1/2. 
— 

Discussion 

By comparing Fig. 2 with Fig. 3, it is apparent that the radial thrust 
is much less efficient than the circumferential thrust for take-off from the 
satellite orbit. For large thrusts, the value of log (% / M; ) for radial thrust 
is more than twice that for circumferential thrust. Furthermore, in case 
of radial thrust, the ratio of thrust to the instantaneous mass, if maintained 
constant, must be larger than As, In case of circumferential thrust, no 
such limit exist. Therefore circumferential thrust is definitely preferred. 

The quantity ሶ log /M,/M,) is a measure of the performance or the 
capability of the vehicle. It has the dimension of a velocity and is actually 
the increase of velocity which the vehicle is capable of in space without 
gravitation. This quantity is conveniently called the characteristic velocity 
of the vehicle. Let this be denoted by V . Then for the case of 
circumferential thrust, Eq. (28) gives 
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where ki is the "escape velocity" from the surface of the earth and À is 


the ratio of the radii of the satellite orbit and the earth. 5 is equal to 


A 1192 km/sec. Fig. 3 then shows that by decreasing the acceleration fe 


EA 


1/2 to 1/3000 g, / , hence the required characteristic velocity V , will 


increase by a factor of two. This is a very important point for the 


designers of space ships. 
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